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Abstract 

^ The dynamical Casimir effect predicts that vacuum amplification effects, resulting in 

the creation of real particles out of vacuum fluctuations, are induced by rapidly mod- 
ulating the boundary conditions of a quantum field. Here we show that a spontaneous 

OO 

release of virtual photons from the quantum vacuum can occur in a quantum optical 
system in the ultrastrong coupling regime. In contrast to the dynamical Casimir ef- 
fect and other pair creation mechanisms, this phenomenon does not require external 
forces or time dependent parameters. It resembles the production of particles during 

O" 1 

the early universe expansion induced by the decay of a false vacuum according to 
inflationary cosmology 
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The dynamical Casimir effect [3] has been recently experimentally realized in a super- 
conducting circuit by modulating the inductance of a quantum interference device at high 
frequencies [3]. Other proposed vacuum amplification mechanisms [I], as the Schwinger 
process [12] and the Hawking radiation [T3], require the presence of huge external fields or, 
as the Unruh effect, the presence of a rapidly nonuniformly accelerating observer [14] . and 
have yet to be observed. The here described spontaneous release of virtual photon pairs in 
the absence of any external drive or observer acceleration resembles the particle production 
during the universe expansion predicted by modern inflationary cosmology. According to 
this theory, a scalar field starting in the false vacuum state eventually decays, and the energy 
that had been locked in it is released to form a hot, uniform soup of particles, which is the 
assumed starting point of the traditional big bang theory [TJ] . 

The Hamiltonian of a realistic atom-cavity system contains so-called counter-rotating 
terms allowing the simultaneous creation or annihilation of an excitation in both atom and 
cavity mode. These terms can be safely neglected for small coupling rates (rotating- 
wave approximation). However, when becomes comparable to the cavity resonance 
frequency of the emitter or the resonance frequency of the cavity mode, the counter-rotating 
terms are expected to manifest, giving rise to exciting effects in cavity QED [61 [151 HE]- 
This ultrastrong-coupling regime is difficult to reach in quantum-optical cavity QED, but 
was recently realized in a variety of solid-state quantum systems joTJTO]. Such regime is 
challenging from a theoretical point of view as the total number of excitations in the cavity- 
emitter system is not preserved, even though its parity is [16j . 

It has been shown that, in the ultrastrong coupling regime, the quantum optical mas- 
ter equation fails to provide the correct description of the system's interaction with reser- 
voirs [IT]. Moreover quantum optical normal order correlation functions fail to describe 
photodetection experiments for such systems [18j [19]. Specifically, for a single mode res- 
onator, the photon rate that can be detected by a photoabsorber is no more proportional 
to (as(t)a(t)) (where a and are the photon destruction and creation operators) but to 
(X~ (t)X + (t)), where X + (t) is the positive frequency component of the quadrature oper- 
ator X(t) = a(t) + a)(t) [19]. The most puzzling property of these systems is that their 
ground state is a squeezed vacuum containing correlated pairs of cavity photons [20]. The 
photon pairs in the ground state |0) are, however, virtual and cannot be detected [21], being 
(0|X~(t)X + (t)|0) = [19J. Otherwise, an observation of a stream of photons from such 
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a system in its ground state would give rise to perpetuum mobile behaviors. Nevertheless, 
we show that such virtual photon pairs, dragged by spontaneous decay processes reach- 
ing or starting from states of the quantum emitter not coupled with the resonator, may 
be spontaneously released and directly observed. A detailed analysis, carried out below, 
shows that such puzzling phenomenon is traceable to a general feature of open quantum 
systems: spontaneous transitions induced by a reservoir occur among eigenstates of the 
total Hamiltonian (environment induced superselection of energy eigenstates [22J ) . Such 
general feature, in presence of a ground state with virtual excitations (or also excited states 
that are not orthogonal to the zero-excitation state), induces the spontaneous Casimir ef- 
fect here presented. A mechanism for the generation of quantum vacuum radiation, based 
on the presence of counter-rotating terms in the light-matter interaction Hamiltonian, has 
already been proposed [151 UHl E3] . Nevertheless these vacuum amplification proposals, as 
the dynamical Casimir effect, require a fast time-modulation or the sudden switch on or off 
of the vacuum Rabi frequency. 

The most promising candidates for an experimental realization of the proposed setting 
are superconducting quantum circuits |24j and intersubband quantum well polaritons [5]. In 
particular, ultrastrong coupling and cascade multi-level structure (the two required features) 
have been experimentally demonstrated in separate experiments with quantum circuits [6j 
[2H 126] and in a single device with multiple quantum wells [5]. 

In the absence of losses, the quantum system depicted in Fig. la,b is described by the 
total Hamiltonian H = H s + i?R a bi with 

#Rabi = U} Q a)a + uJ a cr aa + fi R (a + a t )((T cg + cr ge ) , (1) 

a=g,e 

and H s = u s a ss , where Uq^) (a = s, g, e) are the bare energies of the cavity mode and of 
the atomic-like levels, and o a p describes the transition operators involving the levels of the 
quantum emitter. It is useful to label the eigenstates \j) with increasing eigenenergies flj 
(j integer). Actually the Hamiltonian is block-diagonal and its eigenstates can be separated 
into (i) a non-interacting sector |s, n) with energy u s + nuo, where n labels the cavity photon 
number; and (ii) the dressed atom-cavity states resulting from the diagonalization of -£/R a bi : 
\j) with energy wj. The ground state of this sector can be expanded as 

oo 

|0> = 5>!Uls> 2/c > + <2 fc+ i|e, 2* + 1» . (2) 

k=0 
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FIG. 1: (color online) Sketch of the system under consideration. Two-levels of a single 
quantum emitter are ultrastrongly coupled to a single cavity mode. In the upper panels spontaneous 
emission processes for two different configurations are sketched neglecting the ultrastrong coupling. 
Panels c and d show the involved total energy levels as a function of the coupling Hr/wo and a 
sketch of the allowed spontaneous transitions in the ultrastrong coupling regime, corresponding 
to the configurations in panels a and b respectively. The bending lines as a function of S7r/o;o 
describes the dressed energy levels to--. 

For Q r /lj <C 1, the Jaynes-Cummings (JC) model is recovered and only c° 7^ 0. For 
increasing couplings, e.g. Q-r/uo ~ 0.1, the one photon-pair amplitude c° 2 ~ ^r/^o becomes 
non-negligible. Analogously the excited states can be expanded as 

00 

G>=£(£js.*>+<£*M>)- (3) 

k=0 

In contrast to JC model, they (only those with the same parity of the state \g, 0)) are not 
orthogonal to |g, 0) (being c^ 7^ 0), hence they have a nonzero probability to have zero 
bare excitations. All subsequent calculations are performed at zero detuning: oj — u eg = 0, 
where Uij = Ui — Uj . 

In the first of the two configurations (see Fig. la,c) that we consider, the system is ini- 
tially prepared, e.g. by adiabatic excitation, in the lowest energy dressed state |0) of the 
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system. Although in the ultrastrong coupling regime this state contains virtual cavity pho- 
tons «6[a+a|0> ^ 0), it is a vacuum state for the physical cavity photons that can be detected: 
X + |0) = and hence (0|X~X + |0) = 0. Since |0) is not the lowest energy vacuum, it can be 
considered as the photonic false vacuum, since in this case the state |s,n = 0) is the lowest 
energy (true) vacuum state. Interesting theoretical studies of quantum dynamics in these 
cavity QED systems with |25j and without [27] the rotating-wave approximation recently 
appeared. Panel lc displays the lowest energy levels (the eigenvalues of H: u s + nu and 
ujj) as a function of Or. Let us discuss the spontaneous decay of the initial state |I a ) = |0), 
keeping in mind that spontaneous transitions induced by a reservoir occur among eigen- 
states of the total Hamiltonian (environment induced superselection of energy eigenstates 
[22]). For zero or small coupling rates Or, the initial state |0) reduces to |g, 0) and standard 
spontaneous emission of a photon (at energy u gs ) in the external electromagnetic modes, 
associated with the emitter transition |g) — > cr sg |g) = |s), occurs at a rate 7 gs (fixed by the 
dipole moment of the transition). In the ultrastrong coupling regime, the initial state is 
|I a ) = |0) and possible final states are |Fj) = |s, 2i), where (Fj|cr sg |I a ) = c® 2i . For i = 0, the 
final state contains no cavity photons as in ordinary spontaneous emission. For coupling 
rates Qr/uq < 1, the contribution with i = 1 provides the next dominant term. Hence the 
spontaneous emission of a photon not in the cavity mode at a rate 7 gs comes together with 
a flux of cavity photon pairs at a rate ~ 7g S |c g2 | 2 . 

In the second configuration, shown in Fig. lb, the system is initially prepared into the 
cavity-uncoupled level |Ib) = |s, 0), which is at higher energy with respect to state |e). 
In this case, the initial state can be prepared exciting the matter system from a ground 
state |s'), without involving the emitter states coupled with the resonator. However such 
setting works even in the absence of |s'). We consider the case where spontaneous emission 
induces the following transition: |s) — > |g), as the transition |s) — > |e) is forbidden. For 
small coupling rates Q-r/uq <C 1, standard spontaneous emission of a photon at energy u sg 
directly in the external electromagnetic modes at a rate 7 sg would be observed. In the 
ultrastrong coupling regime possible final states are |Fj) = \j), where (Fj |cr gs | I b ) = c^ . 
The coefficient c^ describes the probability amplitude that the state \j) has zero cavity 
photons. For negligible couplings, only c g0 ^ 0, thus the only possible final state would 
be |0). In this case no cavity photons would be observed, since (0|X~X + |0) = 0. In the 
ultrastrong coupling regime, c^ ^ for even excited states (with the same parity as |g, 0)), 
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hence spontaneous transitions |s, 0) — > |Fj) can occur even for j ^ 0. In particular, the 
next largest coefficients are those for j = 3, 4. These excited dressed states emit physical 
photon pairs ((3|X~X + |3) ^ 0). It is amazing and paradoxical that cavity photon-pairs 
are emitted only because there are excited states of the cavity-matter system which have a 
nonzero probability of having zero cavity photons! 

For describing a realistic system, the cavity and atom dissipation channels need to be 
taken into account. Yet, owing to the very high ratio Q-r/uo, the description offered by 
the standard quantum optical master equation can break down, for example producing 
spurious qubit flipping or photon generation, even at zero temperature [17]. Following Refs. 
|17j . we write the Hamiltonian in a basis formed by the eigenstates of H to describe the 
dissipative processes. We choose a T = temperature environment. Yet generalization to 
T ^ environments is straightforward. We thus arrive at the master equation [T7J 125] . 
pit) = i[p(t),H] + ^2 c £ c p(t), where C c is a Liouvillian superoperator describing the cavity 
(c = 0) and the material system losses c = e — > g, and g — > s (configuration a) or s — > g (b) 
(see Methods for details). 

According to the input-output relations, the emitted photon flux out of the resonator is 
proportional to the mean value of intracavity physical photons: 7o(X~X + ). The results of a 
full numerical demonstration including the cavity and the emitter losses is shown in Fig. 2. 
Figure 2a displays the numerically calculated time evolution of (X~X + ) = Tr[X~X + p(t)] 
for configuration (a) (see Fig. la and c) and for different coupling strengths Qr/ojq = 0.3 (red 
line), 0.4 (blue), 0.6 (black). The signal rapidly grows and reaches a maximum value before 
decaying exponentially due to cavity losses. The signal increases with increasing Qr/uq, as 
a consequence of the buildup of c g 2 - Calculations have been performed at zero detuning and 
by using 7 eg = 70 = 7 gs = 2 x 10~ 2 ^r. Panel 2b shows calculations of (X~X + ) for different 
spontaneous emission decay rates (7 gs /^R = 10~ 2 , 1.5 x 10 _2 ,3 x 10~ 2 ,4 x 10~ 2 ) obtained 
by artificially dropping cavity losses (70 = 0). We also used a coupling rate Qn/u = 0.6 
and 7e g = 7 gs = 2 x 10~ 2 Q R . In the absence of cavity losses, the mean photon number 
reaches a maximum value which does not depend on 7 gs . This result puts forward that the 
phenomenon here investigated is intrinsically different from the dynamical Casimir effect 
where the emitted photon rate depends strongly on the modulation frequency. The inset 
in Fig. 2b displays the maximum cavity mean photon number (X~X + )M ax as a function 
of the energy difference Uq obtained using Qh/uq = 0.8 (other parameters are the same as 
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FIG. 2: (color online) Dynamics of released photons and spontaneous emission spectra. 

a: Time evolution of (X~X + ) for configuration (a) (see Fig. la) and for different coupling strengths 
(see text), b: (X~X + ) for different spontaneous emission decay rates 7 gs obtained by artificially 
dropping cavity losses. The inset in Fig. 2b displays the maximum cavity mean photon number 
(X~X + )Max as a function of the energy difference wg. c: Spectrum S(uj) of spontaneously emitted 
photons, d: Time evolution of {X~ X + ) for the configuration (b) sketched in Fig. lb. 

those used for Fig. 2a). Photon pairs are released only when Uq > u s + 2uq because the 
spontaneous transition |0) — > |s, 2) occurs only when the energy of the final state is below 
that of the initial state (see Fig. lc). Increasing Uq a second rung corresponding to a small 
increase of (X~X + )M ax can be observed. It comes from the activation of the transition 
|0>^|s,4>. 

The release of virtual photon pairs, present in the photonic false vacuum |0), satisfies 
energy conservation. The energy of the radiated photons, induced by the spontaneous decay 
of the false vacuum, comes at the expense of the energy of the quanta emitted into the 
reservoir (e.g. the spontaneously emitted photons in the external electromagnetic modes). 
In the dynamical Casimir effect, the energy of radiated photons comes at the expense of 
the mechanical energy of the moving mirror, which then experiences a friction force by the 
quantum vacuum via the so-called back-reaction effect [271 EE]- In the present case, back- 
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reaction effects can be investigated by calculating the emission spectrum of spontaneously 
emitted photons: S(u) = 1/(2tt) JZc dt IZo ^'(^W^^'))^ 1 ^* , where a ± are the fre- 
quency positive (negative) components of the polarization operator a gs + a sg . The spectrum 
S(oj) of spontaneously emitted photons, calculated for Qr/cjo = 0.6 (other parameters are 
the same as those used for Fig. 2a), is displayed in Fig. 2c. In the absence of interaction 
(dashed line), the spectrum consists of a single Lorentzian peak centered at energy u; gs . In 
the ultrastrong coupling regime the main peak is red-shifted at energy u^ s and a second peak 
at lower energy, centered at Uq s — 2u Q (see Fig. lc) appears. It corresponds to spontaneously 
emitted photons of lower energy associated to the transition |0) — > |s, 2). 




Jet y c t 

FIG. 3: Photon statistics, a: Variance of the squeezed quadrature X = a + a) for different 
coupling rates (see text), b: Normalized second order correlation function g( 2 \t). c: Q^(t,t) = 
(X~ (t)X + (t)) g( 2 '(t,t). d: Normalized third-order correlation function g( 3 \t). 

Numerical results on the second configuration (sketched in Fig. lb) are shown in Fig. 2d. 
The time behaviour of the released photon flux is analogous to that obtained in Fig. 2a for 
the first configuration. The red line has been obtained for u s = 2uq when only the transition 
|s,0) — > |3) is activated (u>^ < u s < u>i). The black line has been obtained for u s = 2.5a>o 
where also the transition |s, 0) — > |4) can contribute. We used ^r/wo = 0.2 and the other 
parameters are the same as those adopted for Fig. 2a. 
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Conspicuous information on the ongoing physics can be obtained studying the statis- 
tics of the emitted photons. A numerical calculation of the variance of the quadra- 
ture X = a + a\ for different coupling rates (Q R /u = 0.3, 0.4, 0.6), is shown in 
Fig. 3a (other parameters are the same as those adopted for Fig. 2a). Squeezing at 
initial times coincides with the squeezing exhibited by the false vacuum state |0) and 
increases with increasing Qr/uo- Nevertheless in the present case squeezing is associ- 
ated to observable photons. The equal-time normalized second order correlation func- 
tion for cavity photons, g^{t) = (X"(t)X-(t)X + (t)X + (t))/(X-(t)X + (t)) 2 , shown in Fig. 
3b certifies a highly super-Poissonian statistics, evidencing that cavity photons are re- 
leased in pairs. Such pair correlation can be further confirmed by calculating Q^ 2 \t,t) = 
(X~ (t)X~ (t)X + (t)X + (t)) I (X~ (t)X + (t)). Such correlation function compares the coinci- 
dence rate with the ordinary photodetection rate. When photons are emitted in pairs 
Q( 2 \t,t) m 1. The decay of these correlation functions is due to cavity losses which tend to 
destroy correlations. Finally we calculated the equal-time third order normalized correla- 
tion function, related to the probability of detecting two cavity photons and a spontaneous 
emitted photon in the external light modes, 

_(8) (f) = (<r-(t)X-(t)X-(t)X+(t)X+(t)a+(t)) 

9 U (a-(t)a+(t))(X-(t)X+(t)) 2 ■ 1 ] 

The obtained huge value of (t, t) (see Fig. 3d) confirms the emission mechanism outlined 
above. Figs. 3b-d have been obtained by using Q^/uo = 0.6 and the other parameters as in 
Fig. 2a. 

In conclusion we proposed a mechanism enabling the release of virtual photon pairs 
which does not require external forces or time dependent parameters. It resembles the 
production of particles during the early universe expansion induced by the decay of a false 
vacuum according to inflationary cosmology. An observation of this mechanism would be 
possible in state of the art experiments with circuit QED and intersubband polaritons. The 
spontaneous Casimir effect opens the way to direct investigation of the most interesting 
feature of ultrastrong coupling, namely the cavity quantum electrodynamic ground state. 
Such vacuum spectroscopy is able to measure the weight of virtual photon pairs in a dressed 
vacuum or the weight of zero photon states in the excited dressed states. Generalizations 
of our study to arrays of coupled cavities would form interesting perspectives for future 
research on virtual photons and pair creation in extended systems (291 E0] ■ In such systems 
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the analogy with the scalar field of the inflationary model would be even more stringent. 
METHODS 

The Liouvillian superoperator in the master equation adopted for the numerical calculations, describing 
the cavity (c = 0) and the material system losses c = e — > g, and g — > s (configuration a) or s — > g (b), is 
defined as C cP (t) = J2j,k>j r i fc, °[|j><*IM*)> with the dissipator V[0]p = \{20pO^ - pO^O - O^Op). Stan- 
dard dissipators are recovered in the limit g — > 0. The relaxation coefficients T 3 c k — 2-Kd c {Akj)a 2 c {Akj)\C^ k \ 2 
depend on the spectral density of the baths d c (Akj) and the system-bath coupling strength a c {A k j) at the 
respective transition frequency A k j = f2fc — £lj as well as on the transition coefficients Cjk = (j\(c+ & )\k) 
[T71HS]. These relaxation coefficients can be interpreted as the full width at half maximum of each \k) — > \j) 
transition. For the sake of simplicity we assume spectral densities d c {A k j) and coupling strengths a 2 to be 
constant. Hence the relaxation coefficients reduce to T^ fe = 7 c |Cj : / J 2 , where j c are the standard damping 
rates of a weak coupling scenario. 
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